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Abstract 

We study the simulation of a single qubit rotation and Controlled-Not gate in a solid 
state one-dimensional chain of nuclear spins system interacting weakly through an Ising 
type of interaction with a modular component of the magnetic field in the z-direction, 
characterized by B z (z,t) = Bo(z) cos St. These qubits are subjected to electromagnetic 
pulses which determine the transition in the one or two qubits system. We use the fidelity 
parameter to determine the performance of the Not (N) gate and Controlled-Not (CNOT) 
gate as a function of the frequency parameter 5. We found that for |5| < 10~ 3 MHz, these 
gates still have good fidelity. 
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1 Introduction 



Almost any quantum system with at least two quantum levels may be used, in principle, 
for quantum computation. This one uses qubits (quantum bits) instead of bits to process 
information. A qubit is the superposition of any two levels of the system, called |0) and |1) 
states, \& = Cb|0) + C\\\) with | C 1 2 + \ C\\ 2 = 1. The tensorial product of L-qubits makes up 
a register of length L, say \x) = ■■■,io), with ij = 0, 1, and a quantum computer with 

L-qubits works in a 2 L dimensional Hilbert space, where an element of this space is of the 
form ^ = ^2C x \x), with ^ \ C X \ 2 = 1. Any operation with registers is done through a unitary 
transformation which defines a quantum gate, and one of the must important result about 
quantum gates and quantum logical operation is that any quantum computation can be done 
in terms of a single qubit unitary operation and a Controlled-Not (CNOT) gate Although 
quantum computers of few qubits [3]- [9] have been in to make serious computer calculations 
one may requires a quantum computer with at least of 100-qubits registers, and hopefully this 
will be achieved in a future not so far away. One solid state quantum computer model that 
has been explored for physical realization and which allows to make analytical and numerical 
studies of quantum gates and protocols [10] is the one made of one-dimensional chain of nuclear 
spins systems 

|11]-|12] inside a strong magnetic field in the z-direction (with very strong gradient in that 
direction) and an RF-field in the transverse direction. Such a model physically is unlikely 
to be constructed, however this represents a good approximation for simulation of quantum 
algorithms and gates whose respective results could be applied in more realistic quantum com- 
puters. Furthermore, the approach relies in the universal character of Quantum Mechanics. 
In this model, the Ising interaction is considered among first and second neighbor spins which 
allows to implement ideally this type of computer up to 1000-qubits or more |13|.|14j. Among 
other gates and algorithms [15], one qubit rotation and CNOT gates were study with this 
quantum computer model [16] . One of the important statement of this model is that one 
keep constant the magnetic field in the z-direction at the location of each qubit. However, 
this statement may be not so realistic in practice for this model or other solid state quantum 
computer based on spin system with very strong axial magnetic field, and then we wonder: if 
there is a magnetic field modulation where this field change slowly with time, how these basic 
elements, one qubit rotation and CNOT gates, would be Of course, in this case, the usual 
analytical approximation without field modulation is not valid anymore, and a full numerical 
calculation is required to see the possible effect of this modulation on 1-qubit rotations and 
CNOT gates. 

In this paper, we want to study this modulation effect of the magnetic field on the Not 
(particular case of 1-qubit rotation, or unitary operation) and CNOT quantum gates. To do 
this, we will assume an additional cosine time dependence on the normal z-direction of the 
magnetic field and will determine, using the fidelity parameter, the minimum variation in the 
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frequency of this modulation to keep these quantum gates elements still well defined. 



2 Quantum Not-gate 

Consider a single paramagnetic particle with spin one-half in a magnetic field given by 

B = (B a cos(cjt), — B a sin(u;t), Bq(z) cos St) (1) 

where the first two components represent the RF-field, and the third component represents the 
strong magnetic field in this direction. The interaction between this particle and the magnetic 
field is given by the Hamiltonian H = —fi- B, where \x is the magnetic moment of the particle 
which is related with the nuclear spin S = HI as = jhl, with 7 the gyromagnetic ratio of 
the particle. So, the Hamiltonian is 

H = —jl- B = -hw Q cos St I z - ^(I + e iu,t + Le^*) , (2) 

where oj q = ^Bq(z ) {z q is the location of the particle) is the Larmor frequency, Q, = ^B a is 
the Rabi frequency, and I± represents the ascent (descent) operator, I± = I x ± il y . If |0) and 
|1) are the two states of the spin one-half, one has that 

Iz\i)= ( ^\i), U\0) = \1), /-|1> = |0). (3) 



To solve the Schrodinger equation, 



ih?^± = H\y), (4) 



one proposes a solution of the form 

|*) = c(<)|0) + Cl (t)|l) (5) 

such that |c Q | 2 + |ci| 2 = 1 at any time. Doing this, one gets the following ordinary differential 
equations 

bJo COS St Q ,-, , . , . 

ic = °— c - - Cl e lu>t (6a) 

and 

uj cos St _ jc . t s 

ici = +^ ci - -c e ^ . (6b) 
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Choosing co(i) = e lujt ^ 2 do(t) and c\(t) = e lult l 2 di(t) in above equations, one has 

9. 



and 



■ oj — uj cos St 
ld o = H 5 d ° ~~ "2 



■ Ct? — w cos St £1 

idi = - di - —do 



(7a) 



(7b) 



which, in turns, can be written as the following uncoupled similar Mathieu equation |19j . 

d + a(t)d = (8a) 
where the complex function a(t) is given by 



«(*) = ^ 

and d\ is obtained from (7a), 



0J n 



ir + oj z 1 — - cos fft 



+ z sin or , 

2 



OJ — UJ cos 

n 



d 



(8b) 



(9) 



For 5 = and on resonance (w = o; ), one has that a = f2 2 /4 , and the system oscillates be- 
tween the states |0) and |1) with and angular frequency corresponding to the Rabi frequency 
fi, as one expected |16| . For 5^0 the solution of this equation is far to be trivial, and instead 
of solving the Eq. (8a), we will find directly the numerical solution of the system (7) with 
the given initial conditions. By taking uj = uj q (resonant case), one expects to obtain the 
transition |0) < — > |1) and to get the quantum Not-gate with a phase. 



To study the performance of the quantum Not-gate as a function of the modulation frequency 
5 , we will calculate the fidelity parameter at the end of a 7r-pulse and make the comparison 
of the ideal wave function, ^expected, with the wave function resulting from our simulation, 

^ sim- 

F = i^sim^ expected) , (10) 

where \^ s im) is the state obtained from numerical simulations, and \^ expected) is the ideal 
expected state, for the initial condition |^ ) = |0), of course, the fidelity coincide with the 
coefficient |ci| 2 . At this point we want to stress that we define \F\ 2 in this way due that 
any quantum gate or algorithm is represented by the final wave function of the quantum 
system. Ideally, if the quantum gate is fully realizable this wave function is represented by 
I ^expected)- However, the non resonant transitions and the error systems (modulation) make 
that the resulting wave function of the complete simulation is given by \^ s i m ). In this way, 
the fidelity is a measure of the good operation of gates and algorithms. On the other hand, 
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there is another measurement for the calculation of the the distance between two states and 
this is the so called Uhlmann- Josza fidelity [T7] . However, in Ref. [18] it has been shown that 
Eq. (10) is a lower bound for the Uhlmann- Josza fidelity. Such a result favors the present 
results. 

Fig. la and Fig. lb show the behavior of the fidelity and the probabilities as a function 
of the parameter 5 at the end of a 7r-pulse, r = n/Cl. We have used the parameters (units 
2ir MHz) = 0.1 and u = 200. The RF- frequency has been chosen equal to the resonant 
frequency u = uj . As one can see, for 5 < 0.2 x 10~ s MHz we can have a very well defined 
quantum Not-gate. 



3 Two qubits model and quantum CNOT gate 

Fig. 2 shows two paramagnetic nuclear particles of spin one-half (qubits) subjected to a mag- 
netic field of Eq. (1), making and angle cos 6 = V3/2 to eliminate the dipole-dipole interaction 
between them. The interaction of the magnetic field with the qubits is carried out through the 
coupling with their dipole magnetic moment /2j = jSi (i = 1, 2), where 7 is the gyromagnetic 
ratio and Sj is the spin of the ith-nucleon (S = HI). The interaction energy is given by 



h = -jt 1 -B 1 -jh-B 2 + hJijpijj s > 



where J is the coupling constant of interaction between nearest neighboring spins, $7 = jB a is 
the Rabi frequency, Hq is the part of Hamiltonian which is diagonal in the basis {|v o )}i =0,i 
and is given by 

H = -h fwi/W +uj 2 ii 2) ) cos eft + HJlPlj®. (12) 



where Wj are the Larmor's frequencies which are defined as 

0Ji = jB ( Zi ) * = 1,2 (13) 

with Zi being the z-location of the ith-qubit. The eigenvalues of Hq on the above basis for 
5 = are 

11 11 

-^00 = +U2- E 01 = -Tri^l - ^2 + 7; J) 

1111 K ' 

By doing co = (En — E\q)/H = 002 — J/2, one gets the resonant transition which defines 
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the CNOT operation 1 10) < — > |11) with a phase involved (e" 1 ^ 2 ), where the left qubits is the 
control and the right one is the target. To solve the Schrodinger equation, 

iH^ = H\V), (15) 
we can assume that the wave function can be written as 

* = C 00 (i)|00) + Coi(t)|01> + Cio(t)|10) + Cu(t)|ll) (16) 
such that J2\ C ij\ 2 = L Thus > 

we arrive to the following system of complex-couple ordinary 

differential equations 

iCoo = ~\ (fa + wa) cos St - \ J) C 00 - | (C i + C 10 ) (17a) 

«Cbi = -i (V - w 2 ) cos^t + i Coi - ^ (Cooe"^ + Cue**) (17b) 
<Cio = (fa - u x ) cos St + ijj Cio - ^ (Cfaoe-^ + C lie !Wt ) (17c) 

iCii = -5 f -(wi + wa) cos St - ^ Jj C n - ^ (C„i + C 10 ) e"*". (17d) 

Doing the transformation C 00 = e iujt / 2 D 00 , C 01 = e-^Ooi, C 10 = e^^Dw, and Cn = 
e~ z3ujt / 2 Du, one gets rid of the fast oscillations and gets the following equations for the 
coefficients D's: 

Woo = + u 2 ) cos St - ^J-uj Dqo - 

ii) 01 = ~ (fa - u 2 ) cos St + ij + u?j D i - 
it) 10 = -i (fa -ui) cos St + + D w - 

ib\\ = -i f — (wi + w 2 ) cos St + Dii 

We solve numerically these equation, and for S = and cj = W2 — J/2, a full transition will occur 
between the states 1 10) and |11). Note that one has CV,(0) = Aj(0) and \Cij(t)\ 2 = \Dij(t)\ 2 . 
For S 7^ 0, we consider two initially conditions cases: Digital case, where the initial condition 
is given by 

|*o) = |10) , (19a) 



^ (D01 + D w ) 


(18a) 


^(D 00 + D n ) 


(18b) 




(18c) 


■ " (£>oi + £>io) • 


(18d) 



G 



that is Coo(O) = 0, Coi(O) = 0, Cio(O) = 1, Cn(0) = 0. Superposition case, where the 
initial condition is 



*o} = \/4l 00 ) + 




y/10 



1 



01 ) + \/T7Tl 10 ) + 




i 



ii) . 



(19b) 



For our simulation, we use the following parameters (units 2ir MHz) = 0.1, uj\ = 100, 
il>2 = 110, and J = 10. The RF- frequency chosen is the resonant frequency to = 102 — J/2, an d 
applying a 7r-pulse, r = ir/Q, we should get the respective CNOT transition 1 10) < — > |H)- 
Fig. 3 shows the behavior of the probabilities and the fidelity as a function of the parameter 5 
at the end of the 7r-pulse and for the digital case. Fig. 4 shows the same as before but for the 
superposition case. This case is more stable (the fidelity decays more slowly than the digital 
case) due to non zero contribution to the terms Coo an d C01 which always contribute with the 
same constant probability 3/10. 
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4 Conclusion 



For a quantum computer model of a chain of qubits in a magnetic field where its z-component 
varies with respect the time, we have studied the Not and Controlled-Not gate behavior as a 
function of the frequency 5 of variation of this component. In general, one can say that for 
5 < W~ 3 MHz these quantum gates remain well defined with a fidelity very close to one. This 
small value in <5 means that it is enough to consider a first order in taylor expansion of the 
cosine function in Eq. (1). We have seen that the fidelity for the superposition case is more 
stable than the digital case due to the contribution to the fidelity parameter of the other no 
zero states involved in the dynamics. Of course, this safety region, defined by 6, for these 
quantum gates does not mean safety for a full quantum algorithm, which is under studied. 
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Figure 1: Quantum Not-gate: (a) Global behavior (b) Local behavior with respect to 8. 
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Figure 4: CNOT behavior, superposition case. 
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